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WHEN CAN GRAVITATIONAL- WAVE OBSERVATIONS DISTINGUISH BETWEEN BLACK HOLES AND NEUTRON 
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ABSTRACT 

Gravitational-wave observations of compact binaries have the potential to uncover the distribution of masses 
and spins of black holes and neutron stars in the universe. The binary components' physical parameters can 
be inferred from their effect on the phasing of the gravitational-wave signal, but a partial degeneracy between 
the components' mass ratio and their spins limits our ability to measure the individual component masses. At 
the typical signal amplitudes expected by the Advanced Laser Interferometer Gravitational-wave Observatory 
(signal-to-noise ratios between 10 and 20), we show that it will in many cases be difficult to distinguish whether 
the components are neutron stars or black holes. We identify when the masses of the binary components could 
be unambiguously measured outside the range of current observations: a system with a chirp mass Ai < 0.871 
M would unambiguously contain the smallest-mass neutron star observed, and a system with A4 > 2.786M Q 
must contain a black hole. However, additional information would be needed to distinguish between a binary 
containing two 1.35 M Q neutron stars and an exotic neutron-star-black-hole binary. We also identify those 
configurations that could be unambiguously identified as black-hole binaries, and show how the observation of 
an electromagnetic counterpart to a neutron-star-black-hole binary could be used to constrain the black-hole 
spin. 

Subject headings: black hole physics — gravitational waves — stars: neutron 



1. INTRODUCTION 

The first observing runs of the Advanced Laser Interferom- 
eter Gravitational-wave Observatory (aLIGO) are expected in 
~ 2 015, with Advanced Virgo following on a similar sched- 
ule dHarrvet al.l 120 10t lAcernese et"aT] 120091) . The primary 
source for these observatories is the coalescence of binaries 
containing black holes and/or neutron stars, with predicted 
rates between a fe w and several hundred per year at detecto r 
design sensitivity ( Aba die et al.ll2.Q10l : iDominik et alj|2012l) . 
The components of these binaries are formed in supernovae 
when the core of a massive star collapses to a compact rem- 
nant, although the exact collapse mechanism remains un- 
known. Detailed knowledge of the mass distribution of black 
holes and neutron stars will provide vital clues to their forma- 
tion as well as explore the equation of state (EOS) of nuclear 
matter at high densities. 

Measuring the upper and lower limits of neutron-star 
masses allows us to constrain the super nova engine and the 
nuclear physics of neutr on-star remnants (Lattimer & Prakashl 
120 lOt iFrver et all 120 1 2h . For example, the collapse of low- 
mass stars (8- 12M Q ) is believed to quickly produce explo- 
sions with very little mass accreted in a convective engine 
phase. The predicted mass of the compact remnant will be 
less than the Chandras ekhar mass by an amount that depends 
on the collapse mo del ( IFrver et al.ll 19991: iKitaura et aLll2006t 
lDessartetal.ll2007h . If we can place an upper limit on the 
mass of a low mass neutron star (in the 1.0- 1.2M Q range), 
we can distinguish between current models, effectively us- 
ing these low-mass systems to constrain the EO S and neu- 
trino p hysi cs in core- c ollaps e. At higher masses, lOzeT et al. 
and iFarr et all (1201 11) have argued that there is a gap 
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between ~ 2-4M© where no compact objects exist. If 
true, the amount of material that falls back onto the newly 
formed compact remnant must be small, arguing against en- 
gine mecha nisms that take a long time (more than 200 ms) 
to develop (Belczvn ski et al.|[2012l) . However, this mass gap 
may simply be an artifact of poor mass resolution of X-ray bi- 
naries and poor statistics (Kr eidberg et alj|2012l) . Black-hole 
mass distributions will allow us to explore the fall-back of ma- 
terial in a weak supernova explosion and black-hole masses in 
solar metallicity environments will provide clues into stellar 
mass loss. With an accurate black-hole mass distribution, we 
can study these open questions in stellar evolution. 

The binary's gravitational-wave phasing depends at 
leading order on its chirp mass M. = {m\miy'l i {m\ + 
w^)" 1 / 5 , where m\ and rri2 ar e the binary's component 
masses dPeters & Mathewslll963l) : this quantity will be most 
accurately measured in a gravitational-wave detection. The 
mass ratio r\ = m\m2/(mi+m2) 2 enters through hig her-order 
correc tions and is less accurately measured; see e.g.. lBlanchetl 
( 120021) . There is also a partial degeneracy between the mass 
ratio and the angular momentum \ l .2 = J\ .il m \ 2 of each com- 
pact object (the spin), which further limits our ability to mea- 
sure the binary's component masses. Heuristically, this can 
be understood as follows: a binary with spins aligned with 
the orbital angular momentum will inspiral more slowly than 
a non-spinning system. Similarly, a binary of the same to- 
tal mass but with more extreme mass ratio will inspiral more 
slowly. However, a binary with the same chirp mass but with 
more extreme mass ratio will inspiral more quickly. The ef- 
fect on the waveform of decreasing r} can be mimicked by 
increasing the component spins. 

We investigate the accuracy with which the component 
masses can be determined from gravitational-wave observa- 
tions of binary-neutron-star (BNS), neutron-star-black- hole 
(NSBH), and binary-black-hole (BBH) systems, focussing on 
systems where the object's spins are aligned with the orbital 
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angular momentum. Since the first signals detected by aLIGO 
are likely to have signal-to-noise ratios (SNRs) close to th e 
observable network threshold of ~12 ( lAbadie et alj|2012h . 
we focus on signals with SNRs 10-20, which will account 
for ~80% of observations. For these SNRs, we find that the 
mass-ratio-spin degeneracy will prevent us from accurately 
measuring component masses. We identify the region of the 
mass parameter space for which it will not be possible to de- 
termine whether the compact objects are black holes or neu- 
tron stars using gravitational-wave observations alone, when 
we can conclusively measure compact-object masses outside 
the currently observed limits, and show how the observation 
of an electromagnetic counterpart to an NSBH could be used 
to constrain the black-hole spin. 

2. PARAMETER ESTIMATION METHOD 

Fisher-matrix methods show that the binary's chirp mass 
is recovered well by matched filtering, w ith accuracies of 
~ 0.0 1% for typic alBNS systems in aLIGO dFinn & Chernoffl 
[T99llArun etaDl200lh . If we assume that the neutron stars 
are n onspinning 77 can be measured to an accuracy of <~ 
1.3% dArun et al.ll2005l) . Estimates of the effect of the mass- 
ratio-s pin degeneracy were first m ade by ICutler & Flanagan! 
d 19941) and lPoisson & Willi d 1995b using the Fisher approach. 
The degeneracy between the mass ratio and the total effec- 
tive spin x = ( m iXi + m 2X2)/( m i +mi) degrades the ability 
to measure the mass ratio and hence the component masses. 
We go beyond the se studies, using the method introduced in 
iBaird et all (120121) to equate a confidence interval with a re- 
gion where the match between the signal and model wave- 
forms exceeds a given threshold. We use this method to in- 
vestigate parameter degeneracies for a wide range of binaries 
and interpret the expected measurement accuracy in the con- 
text of the astrophysical questions discussed above. 

We model the w aveforms with the TaylorF2 inspi- 
ral approximant (ISathvaprakash & Dhurandharl |1991[ 
ICutler & Flanagan! 119941: lDrozetal.1 I1999F to leading 
order in amplitude and 3.5 p o st-Newtonian (PN) order 
in phase dBlanchet et alJ [19951 l200l l200l 120051) with 
spin-o rbit terms to 3.5PN order and spin terms to 2PN 
order dKidder et al .11 19931: lKidderf[T995l) . For systems with 
total masses below ~ 8M , our results with TaylorF2 
are consistent with those from phenomenolo gical BBH 
mode l s that include the mer ger and ringdown dAjith et aD 
1201 It ISantamaria et al.ll2010h . calibrated agai nst numerical- 
relativ ity waveforms with mass ratios up to 1 :4 (lHannam et al.l 
120101) . For the higher-mass BBH results in Sec. |4]we use the 
full merger model. Throughout, we assume for simplicity 
that the component spins are aligned with the orbital angular 
momentum. In this case, the binary's distance, orientation 
and sky location affect only the overall amplitude of the 
waveform, and we do not consider them here. 

For two waveforms hi and /12 the match is given by 



M= max 
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where (a\b) is the standard noise- weighted inner product 
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In all cases we use a noise sensitivity S„(f) corre- 
sponding to the zero-detuned high-power configuration of 



aLIGO dShoema ker 2009) with a 15 Hz low frequency cut- 
off. We construct a 90% confidence region for a signal 
in the (mi,ni2,x) space, which corresponds to the three- 
dimensi onal region where M > 0.968 (0.992) for an SNR of 
10 (20) dBaird et al.ll2012l) . This method is more accurate at 
low SNRs than the Fisher-matrix approach. 

3. NEUTRON-STAR BINARIES 

Observed neutron-star masses curr ently lie be- 
tween 1.0 ± O.IOMq dRawlsetal.1 1201 lb and 
1.97±0.04M Q dDemorest et al.1 120 101) . General relativ- 
ity and causality place a strict u pper limit on the m a ximum 
neutron-star mass of 3.2M dRhoades & Ruffinil 119741) : 
the actual maximum mass is determined by the as yet 
unknown neutron-star EOS. The observed masses of stars 
in double neutron-star systems is narrower with a pe ak at 
1.35M and a width of 0.13M Q dKiziltan et alj|2010l) . We 
begin by considering a canonical BNS system with masses 
mi = ni2= 1.35M Q and no spins. Fig. Q] shows the regions 
of the mass plane that are consistent with this source at 
90% confidence for SNRs 10 and 20. The component 
masses consistent with the signal lie roughly along a line of 
constant chirp mass. However, there is significant spread 
in the recovered component masses due to the degeneracy 
between 77 and \ 1.2 in the gravitational-wave phase evolution. 
Typically, the degeneracy persists over a range of 0.3 in x- 

If we assume that the compact objects are non-spinning, the 
component masses are recovered in the ranges 1 .15-1.35M 
and 1.35-1.6M0. The observed spins of double neutron stars 
are low, with t he minimum obser ved period of 22.70ms for 
J0737-3039A dBurgay et al.ll200l . i.e., X ~ 0.02, where the 
neutron-star period can be related to the spin by approxi- 
mately 

-(S0(?H(T)- 

If we constrain the neutron stars to have a spin \ < 0.05, 
the range of consistent component masses extends to 1.0- 
1.35M and 1.35-1.9M . Therefore at typical SNRs, we 
would be unable to distinguish the canonical 1.35M BNS 
from a more exotic BNS with m 1 = l.OM ,ra2 = 1.9M . 

The fastest spinning pulsar (PSR J1748-2446ad ) has a 
period of 1.4ms, or \ ~ 0.3 dHessels et alJ 120061) . Neu- 
tron stars are considered unlikely to have a period less than 
1 ms (jChakrabartv 20081). a lthough breakup frequencies could 
be a factor of - 2 higher (lLo & Linll20Tlb . If we allow for 
larger component spins, the region consistent with the canon- 
ical BNS system extends to 0.65-1. 35M and 1.4-3.1M . 
The maximum total spin of binaries within the 90% confi- 
dence region at SNR 10 is ^0.3; see Fig. Q] inset. Based 
on gravitational wave observations alone, a binary of two 
1.35M neutron stars could not be definitively distinguished 
from a low-mass neutron star and a black hole or neutron star 
in the mass gap, even at SNR 20. If the larger object is a neu- 
tron star its spin would be unusually large, but not impossible. 
To constrain the larger object's mass below 2 M without as- 
sumptions on the spins would require a signal with an SNR of 
>40, expected for ~ 2% of observations. 

In Fig. 12 we consider a family of equal mass, non-spinning 
BNS systems with component masses at the lowest (l.OM ) 
and highest (2.OM ) observed neutron-star masses, and at the 
highest theoretical neutron-star mass (3.2M ). The 1.35M 
binary is shown for reference, and the observed mass gap is 
indicated by the shaded region. Again, the 90% confidence 
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FIG. 1.— The 90% confidence region around a 1.35M -1.35M Q BNS 
system at SNR 10 (gray) and SNR 20 (black) in the zero-detuned high-power 
configuration of aLIGO. The vertical bars indicate where the SNR 10 confi- 
dence region would be truncated if we restrict to non spinning neutron stars, 
and if we restrict to neutron-star spins less than 0.05. The inset shows the total 
effective spin x = ('« l X 1 + m 2X2)/('"l +'"2) of the waveforms, with respect 
to the mass of the larger body. 



regions follow lines of approximately constant chirp mass. 
We can see that for an exceptionally low or high mass binary, 
we will be able to identify at least one of the components as 
extraordinary. The observation of a BNS with a chirp mass 
M. < 0.871 M Q would yield the unambiguous detection of a 
compact object with mass less than 1.0M©. The less mas- 
sive component could have a mass in the range O.5-1.OM , 
but any mass in that range would challenge our current under- 
standing of neutron stars and their formation in supernovae. 
A chirp mass Ai < 1 .045 would yield an unambiguous detec- 
tion of a < 1 .2M Q neutron star. This constraint would rule out 
many modern calculation s of stellar collapse dKitaura et al.l 
l2006tlDessart et ai1l2007f) and provide a strong validation test 
for future calculations. 

Similarly, the observation of a system with M.>\ .741 M Q 
would indicate the detection of a neutron star of mass > 
2.0 Mq provided we can unambiguously identify both com- 
ponents as neutron stars. However, the observation of a BNS 
system with m\=m 2 = 2.0M© is consistent with a binary con- 
taining a neutron star and a black hole in the mass gap or an 
exotic, low-mass neutron star and a black hole. Although any 
configuration along this line is of interest, we cannot strongly 
constrain the component masses with gravitational-wave ob- 
servations at low SNR. In particular, even if we assume that 
neutron stars are nonspinning and that the minimum BH mass 
is 3M Q , we cannot rule out that the possibility that it is a 
NSBH system: our assumption would only remove the por- 
tion of the confidence interval with 2.4M© <m 2 < 3.0M©, 
and the observation would be of either a BNS system with 
well-constrained masses, or an NSBH. 

Finally, we consider a binary with component masses at the 
upper end of the neutron-star mass limit, above which we ex- 
pect the components to be black holes. In this case, we could 
conclusively say that one of the components must be a black 
hole, although the degeneracy limits our ability to draw strong 
conclusions on the component masses. The binary could ei- 
ther be a BBH system with component masses in the mass 
gap, or an NSBH system with masses consistent with previ- 
ously observed compact objects. 

Note that in all the cases we have considered, the equal- 




FlG. 2. — The 90% confidence regions for a number of different non- 
spinning compact-binary configurations; see text for interpretation. The num- 
ber at the end of each confidence region is the chirp mass of the binary. The 
gray shaded region indicates the current observational mass gap. 



mass line provides a hard upper (lower) limit on m\ {mi). 

Fig.|2]shows that it will be difficult to explore the mass gap 
with single gravitational-wave observations at low SNRs. To 
constrain one of the masses to lie within the mass gap, the sys- 
tem would need to be observed with an SNR of greater than 30 
(~ 5% of observations); in the case of the m\ = m 2 = 2.0M Q 
binary, the mass on one of the objects would then be con- 
strained between 2M© and 3.5M©, placing it directly within 
the mass gap. 

4. BINARY BLACK HOLES 

Black hole spins can vary between and 1, with obser- 
vatio ns of X-ray binari es supporting the full range of val- 
ues dZhang et al.l 11997b . As before, the mass-ratio-spin 
degeneracy precludes precise measurement of the compo- 
nent masses. However, if we observed a BBH with m\ = 
m 2 = 36 Mq, corresponding to a binary with both compo- 
nents above 35 M Q — the largest ma ss observed for a stellar- 
mass black hole in an X-ray binary (ISilverman & Filippenkol 
120081) — we could conclude that one of the black holes has a 
mass of 36 M Q or higher, providing the first observational ev- 
idence for a stellar-mass black hole above 35 M©. The same 
is true for any binary with the same or larger chirp mass, 
M > 31.34 M©. 

A BBH will only be unambiguously identified if all bi- 
naries within the 90% confidence interval for the measured 
component masses are also BBH systems. The shaded re- 
gion in Fig. [3] shows the part of parameter space in which 
BBH systems could be identified as such at SNR 10, as- 
suming a maximum neutron-star mass of 3.2M©. The up- 
per curve shows the boundary of this region computed using 
the inspiral only (TaylorF2) model. The lower curve is com- 
puted using a phenomenolog ical inspiral-merger-ringdown 
model ( iSantamaria et al.l 120101) . When the additional infor- 
mation from merger and ringdown is included, a larger region 
of the parameter space can be identified as a BBH system, 
illustrating the impact of merger-ringdown in these systems. 

5. ELECTROMAGNETIC COUNTERPARTS 

Fig. [4] shows the 90% confidence regions for an NSBH, 
where the NS has mass 1 .4 M© and no spin, and the BH has 
m 2 = 10M© with spin \i = 0-7 aligned with the orbital an- 
gular momentum. The uncertainty in both masses is large: 
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FIG. 3. — The shaded region indicates binaries would be unambigu- 
ously identified as BBH systems, assuming a maximum neutron star mass 
of 3.2 Mq (dotted line). The light gray region represents the most conserva- 
tive estimate of this region, using an inspiral-only waveform model, while the 
additional darker region uses a full merger-ringdown model. 
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FIG. 4.— The 90% confidence regions at SNRs 10 and 20 for an NSBH 
system with masses mi = 1.4Mq and 1112 = IOMq, and BH spin \2 = +0.7. A 
GRB observation would allow an estimate of the BH mass of mi = 1 1 ± 3M0 , 
and X2 = 0.7 ± 0.05 (see text). 



this could be an NSBH system containing a BH between 
6M Q and 26 M Q , or an equal-mass BBH system, with masses 
m\ = ni2 = 3.5 Mq. The possible mass ratio of this system ex- 
tends from 1:1 through to 1:25. Even at SNR 20 the mass 
ratio ranges from 1:2.4 up to 1:11. We note that the wave- 
form models used here only include the waveform's dominant 
harmonic and do not include merger and ringdown (the phe- 
nomenological model is calibrated only up to mass ratio 1:4). 
The confidence region may significantly change if higher- 
order or merger-ringdown effects are included, emphasizing 
the need for further waveform development in interpreting 
gravitational-wave observations as well as for detection. 

The observation of a gamma-ray burst or other electromag- 
netic signal would unambiguously identify one of the com- 
ponents as a NS and break the mass-ratio-spin degeneracy. 
If we accept th at the mass of a ne utron star in a binary is 
1 .35 ± O.25M dKiziltan et aUHoTot) . then the BH mass is re- 
stricted to mi = 11 ± 3M Q . The further assumption that the 
neutron-star spin is below 0.4 restricts the possible values of 



the black-hole spin, allowing an estimate of \2 = 0.7 ±0.05. 
This is shown in Fig. [4] (inset), which shows the 90% confi- 
dence region restricted to the xi-mi plane, and the region for 
which the neutron-star mass is between I.IMq and 1.6M Q . 
However, due to beaming only a fraction of gravitational- 
wave observations of BNS and NSBH systems are likely to 
be accompanied by GRB observations (lAbadie et alj|2012ab . 

6. CONCLUSIONS 

Although gravitational-wave observations will accurately 
measure the chirp mass M. of binary mergers, a mass-ratio- 
spin degeneracy prevents the component masses and spins be- 
ing measured accurately at low SNRs. In many cases it will 
be difficult to determine whether the components of the binary 
are neutron stars or black holes. However, we have illustrated 
several cases where significant results can be inferred from 
gravitational-wave observations. 

We illustrate several situations where the binary must con- 
tain (at least one) compact object that is more or less mas- 
sive than anything observed to date; M. < 0.871 M Q indicates 
a neutron star smaller than 1 Mq, M. > 1.741 M© a neutron 
star larger than 2 M Q (if we can independently verify that it is 
a BNS system), and M > 31.34M Q a black hole larger than 
35 M Q . Observations at SNRs higher than 30 will be required 
to clarify the existence of the mass gap. The observation of 
an electromagnetic counterpart will, in certain situations, al- 
low us to identify the system as an NSBH and, if we know the 
distribution of neutron-star masses, measure the black-hole 
mass and spin with high accuracy. 

Our results are qualitatively robust, but the true confi- 
dence intervals may vary; they were genera ted by calculat- 
ing waveform mismatches dBaird et al. 12012b. rather than fu ll 
parameter-estimation m e thods (Ivan derSluvs et ai1l200 8a b; 
IVeitch & Vecchiol 120 10t iFeroz et al.l 12009b . In "some cases 
inclusion of merger-ringdown, higher-order post-Newtonian 
corr ections to spin e ffects, particularly at high mass ra- 
tios (lArun et al.ll2.009t) . and tidal effects will slightly alter the 
size of the error regions. The issues only serve to highlight 
the urgency of improved waveform models. The inclusion of 
precession and higher harmonics will modify our results, but 
we do not e xpect these to sig nificantly break the degeneracy 
we discuss (Baird et al. 2012). 

Our results highlight the significance of the mass-ratio-spin 
degeneracy in gravitational-wave observations of compact bi- 
naries. This needs to be explored further using more com- 
plete parameter-estimation methods, more accurate waveform 
models, and extended to include an exhaustive study of the 
parameter space. Results from a population of signals, rather 
than individual observations, should also be investigated. Un- 
derstanding the uncertainties in masses and spin measure- 
ments will be essential to interpreting gravitational-wave ob- 
servations. 
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